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We study the possibility of spin injection from Fe into Si(OOl), using the Schottky barrier at 
the Fe/Si contact as tunneling barrier. Our calculations are based on density-functional theory for 
the description of the electronic structure and on a Landauer-Biittiker approach for the current. 
The current-carrying states correspond to the six conduction band minima of Si, which, when 
projected on the (001) surface Brillouin zone (SBZ), form five conductance hot spots: one at the 
SBZ center and four symmetric satellites. The satellites yield a current polarization of about 50%, 
while the SBZ center can, under very low gate voltage, yield up to almost 100%, showing a zero- 
gate anomaly. This extremely high polarization is traced back to the symmetry mismatch of the 
minority-spin Fe wavefunctions to the conduction band wavefunctions of Si at the SBZ center. 
The tunneling current is determined by the complex band structure of Si in the [001] direction, 
which shows qualitative differences compared to that of direct-gap semiconductors. Depending on 
the Fermi level position and Schottky barrier thickness, the complex band structure can cause the 
contribution of the satellites to be orders of magnitude higher or lower than the central contribution. 
Thus, by appropriate tuning of the interface properties, there is a possibility to cut off the satellite 
contribution and to reach high injection efficiency. Also, we find that a moderate strain of 0.5% 
along the [001] direction is sufficient to lift the degeneracy of the pockets so that only states at the 
zone center can carry current. 

PACS numbers: 72.25.Hg,85.75.-d 



I. INTRODUCTION 

The electrical injection of spin polarized carriers, in 
short called electrical spin injection, into the conduction 
band of semiconductors, is one of the key elements for 
realizing spin transistors such as the one proposed by 
Datta and Dasi Following arguments by Schmidt and 
co-workersj2 Rashba^ and Fert and Jaffres^ it was re- 
alized that efficient spin injection requires the presence 
of a tunnel barrier at the ferromagnct /semiconductor in- 
terface. Significant experimental success started in 2001- 
2002^ in junctions that contained a tunnel barrier be- 
tween the ferromagnetic metal and the semiconductor, 
either in the form of the Schottky barrier of the semi- 
conductor itself or in the form of an ultra-thin film of 
some other insulating material^ such as AlO or MgOii 
Most works have focussed on injection into GaAs^ii 
where the current polarization can be optically detected 
in GaAs/AlGaAs/GaAs quantum wells. Recently, how- 
ever, electrical spin injection into Si has also been demon- 
strated both via electrical spin detection^u^ and optical 
detection^ with impressively large spin coherence length, 
reaching up to 350 /im^ 

The present work was motivated by the increasing ex- 
perimental activity in the direction of spin injection and 
manipulation in Si^^ i 10 ' 11 ^ 2 ' 13 as well as by arguments 
for advantages of spin transport in SiJ^ Based on density- 
functional calculations we examine the possibility of di- 
rect spin injection from Fe into Si, with the Schottky 
barrier of Si used as the necessary tunneling barrier. We 
focus on idealized epitaxial, atomically flat Fe/Si(001) in- 
terfaces. After briefly presenting the magnetic structure 
of the Fe/Si interface, we discuss the tunneling proper- 



ties of Si based on the concept of complex band structure. 
We then show that spin injection is possible and can be 
tuned by adjusting the Fermi level Ep in the band gap 
of Si in the tunneling region, or by straining Si in the 
injection region in order to cut off satellite contributions. 
We further discuss the so-called f -point rule for increased 
injection efficiency and show that, in the case of Fe/Si, 
it is harder to satisfy than in Fe/GaAs. 



II. METHOD AND DETAILS OF 
CALCULATION. LIMITATIONS OF APPROACH 



The electronic structure of the junctions was cal- 
culated within the local (spin) density approximation, 
L(S)DA, to density functional theory The Kohn- 
Sham equations were solved using the screened Korringa- 
Kohn-Rostokcr (KKR) Green function method^ A non- 
relativistic treatment was chosen, as it is known that rel- 
ativistic corrections result in a reduction of the calculated 
semiconductor gaps within the LDA. The atomic sphere 
approximation was used in most cases, except in par- 
ticular tests which were performed with a full potential 
treatment. The conductance was calculated in the zero- 
bias limit within a Landauer-Biittiker approach adjusted 
to the KKR Green function method*^ 

The self-consistent electronic structure of the junction 
was calculated using the decimation technique^ for the 
consideration of the half-infinite leads. An angular mo- 
mentum cutoff of ^ m ax = 2 was used. The details of the 
electronic structure that are important here change only 
little between f max = 2 and € max = 3; for example, the 
spin moment of bec Fe changes by less than 5%. However, 
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FIG. 1: Schematic representation of the band diagram at the 
Fe/Si interface used in the calculations. 



an angular momentum cutoff of t max = 3 was necessary^ 
for the calculation of the conductance matrix elements. 

The Fe/Si(001) interface was taken to be ideal (atom- 
ically flat). Contrary to the wavefunctions, the pertur- 
bation of the potentials dies out relatively fast with the 
distance from the interface. Within a finite region of 
two atomic layers around the interface, the potentials 
are perturbed from their bulk-like values and calculated 
sclf-consistcntly by sandwiching 5 layers of Si between 
two Fe leads. This way also provides the band alignment 
between Fe and Si. Subsequently the half-infinite Si lead 
is constructed by maintaining the potentials only up to 
the middle of the junction (left Fe lead and Si up to third 
layer), and repeating the potential of the third Si layer 
(which is considered already bulk-like) ad infinitum by 
use of the decimation technique. The decimation tech- 
nique is also used to produce the half-infinite, bulk-like 
Fe lead. The same method has been applied in previous 
works on Fe/GaAs and Fe/ZnSe-i^^. 

We consider the following external parameters that 
are involved in an injection experiment (see also Fig. [l}. 
Firstly, the Schottky barrier has a thickness of Nsb layers 
and a height determined by the band alignment of Si with 
respect to the Fermi level of Fe. Experimentally, this can 
be adjusted by appropriate doping or by a gate voltage. 
In the calculations, this is simulated by adding an appro- 
priate constant shift to the Si potentials up to Nq,b, so 
that, within the barrier, the bands arc shifted and Ep 
(which is fixed by the Fe lead) falls at a desired position 
between the valence band and conduction band edges, 
E v and E c : Ep = E v + E tunn . In this step the first two 
Si layers at the interface are excluded, as it is considered 
that their electronic structure will be mainly affected by 
their proximity to Fe. Secondly, after the barrier region 
the band alignment must be such that the injected elec- 
trons arrive in the Si valence band. Experimentally the 
appropriate band alignment can also be adjusted by a 
gate voltage or by doping. Here, it is again accounted 
for by adding a proper energy constant to the potentials 
of all Si layers from Nsb + 1 and on, so that Ep falls 
slightly in above E c : Ep = E c + S. In the present work 
S = 28 meV (2 mRyd) was used except where otherwise 
indicated. 

A few comments are in order on the chosen band di- 
agram. Firstly, in experiment a Schottky barrier will 
possibly have a different, less abrupt shape. We defer 



this question for Sec. IVI1 because after the discussion 
has revealed the effect of the positioning of Ep in the 
gap. Secondly, the chosen value of 5 = 28 meV can be 
perhaps reasonable for GaAs, but would probably result 
in an unrealistically high carrier concentration in Si due 
to the high Si effective mass. This choice helps in vi- 
sualizing and understanding the current distribution in 
the SBZ. In Sec. [V] we also show calculations with S up 
to 1.4 meV, which in fact reveal a qualitative change of 
injection efficiency. 

It is well known that local density-functional theory (as 
is the LDA or the generalized gradient approximation) 
underestimates the gap in semiconductors and insulators; 
in the present calculations the Si band gap is found to 
be approximately 0.38 eV. Thus, the values of tunneling 
conductance obtained within density-functional calcula- 
tions can only provide qualitative understanding, includ- 
ing trends, but not quantitatively correct results. Qual- 
itatively, the underestimated band gap in the tunneling 
region can be partly compensated by increasing the bar- 
rier thickness Nsb- Here, Nsb was varied up to 60 Si 
layers. 

Especially for the calculation of the bulk Si band struc- 
ture (including the complex band structure), significant 
quantitative improvement can be achieved by using an 
orthogonalizcd plane wave method with pseudopotentials 
fitted to optical excitation experiments. Here we followed 
such a method with the parameters taken from Ref . [20l ] ; 
in the presentation of the results we indicate when this 
method has been used. 

Furthermore we note that an atomically flat interface 
is, at this point, an idealization. In an experiment, most 
likely interface roughness and some form of iron silicide 
will be present at the interface. However, it cannot be 
excluded that good quality, atomically almost flat inter- 
faces can, in principle, be made, as has been achieved 
in Fe/MgO junctions. Moreover, the parameter space 
becomes simply too big if many possibilities of interface 
structure are to be calculated. Since modern experimen- 
tal techniques allow for a detailed imaging of the interface 
structure^ it is possible to relate further calculations to 
such experimental input. In Sec. I VII we point out which 
calculated properties are specific to the idealization of a 
flat interface. 

Finally we comment on the lattice mismatch of Fe and 
Si. The Si lattice constant (5.43 A) is approximately 
twice the one of Fe (2.87 A) with a mismatch of about 5%. 
We adopt the Si lattice constant for the calculation (the 
in-plane unit cell accommodates now two Fe atoms per 
Fe layer), to simulate the situation where a thin film of 
Fe is in contact with a thick Si barrier. The stress on the 
Fe contacts will be appreciable, but it is sufficient that Fe 
retains the prescribed epitaxial structure for only a few 
layers for the our conclusions to be correct. This has been 
demonstrated in calculations on Fe/MgO/Fe junctions,— 
where it was shown that the transport properties depend 
on the structure of the Fe lead only close to the interface: 
even if Fe becomes amorphous a few layers away from 
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FIG. 2: (color online) Atom-resolved spin moments close to 
the interface in a Fe/Si/Fe(001) junction, calculated at the 
Si lattice parameter. Evidently the Fe moments at the inter- 
face are strongly suppressed. The two values for each layer 
correspond to the two inequivalent atoms. 



the interface, the spin filtering properties practically do 
not change. However, the choice of lattice constant does 
affect the moments significantly, and therefore tests on 
interface relaxation and tetragonalization were made, as 
we discuss in Sec. IIIIl 



III. SPIN MOMENTS AT THE INTERFACE. 
EFFECT OF LATTICE CONSTANT 

The spin moment of bulk Fe was found not to change 
too much with lattice constant: a moment reduction of 
12% was calculated at the Si lattice constant compared to 
the result at the Fe lattice constant. However, the change 
of the moments at the interface is appreciable. At the Si 
lattice constant, the moments of the interface Fe atoms 
are reduced to approximately half their bulk value, as is 
shown in Fig. [2j This effect, previously reported, e.g., in 
Ref. [lH, is riot present in calculated junctions at the Fe 
lattice constant, where the interface Fe moments are not 
suppressed at all. 

The magnetic structure was further tested by calcu- 
lating a finite slab of 4 layers Fe/ 4 layers Si / 4 layers 
Fe with a full potential treatment and £ max = 3. Here 
the in-plane lattice constant of Si was adopted, while the 
Fe layers were allowed to assume an optimal c/a ratio, 
and the interlayer distance at the interface was relaxed. 
No qualitative differences were found compared to the re- 
sults using the Si lattice constant for a and c: the strong 
suppression of the moments was still present. 

We conclude that a strong suppression of the Fe in- 
terface moment is unavoidable if the in-plane lattice pa- 
rameter of Si(001) is adopted. Interestingly, this does 
not qualitatively affect the transport properties that in- 
terest us here, in particular the spin polarization of the 



injected current. This was found by test calculations of 
spin injection in a Fe/Si(001) junction using the Fe lattice 
parameter (whence the moments at the interface where 
not compromised). 



IV. ROLE OF SILICON REAL AND COMPLEX 
BAND STRUCTURE 

Silicon has an indirect band gap. The valence band 
has a maximum at the T point (k = 0), while the conduc- 
tion band has six equivalent degenerate minima at k = 
(±fc o ,0,0), (0,±fc o ,0), (0,0,±fe o ) with k w 0.85 2n/a 
(a is the lattice parameter). These properties are most 
important for a qualitative discussion of the tunneling 
conductance and electron injection. 

The results of this section were obtained by using 
pscudopotcntials fitted to optical transitions^ The band 
structure of Si around the gap is shown in Fig. [3l The 
conduction band around the minima forms six ellipsoidal 
"pockets" of highly anisotropic effective mass, as shown 
schematically in Fig. [4] Electron injection takes place 
into these pockets; if the junction growth is along the 
[001] direction, then the pockets are projected in five 
conductance hot spots on the two-dimensional (001) sur- 
face Brillouin zone (Fig. 0J. One of these is at fen = 0, 
while the other four are equivalent and form a satellite 
structure. (We denote by fe|| the projection of the Bloch 
fc-vector on the surface Brillouin zone.) 

The pocket structure is highly important also for the 
tunneling properties. Setting the z axis along the [001] 
epitaxy direction, the tunneling wavefunction at Ef de- 
pends on the decay parameter K(kt\;Ep) as ip(r) ~ 

cxp(— kz). The set of n{k\\\E) forms the so-called com- 
plex band structure; evidently, the lowest values of k are 
responsible for the highest tunneling current. The value 
of /c(fc|| ; Ep) depends on the proximity of Ep to the Bloch 
eigenvalues E(k«; k z ) of the same kn in the conduction or 

valence band, and on the curvature of E(ku ; k z ) with re- 
spect to k z (i.e., the effective mass in the k z direction). 
Thus the effective mass anisotropy of the pockets in the 
Si tunneling barrier has consequences that we now dis- 
cuss. 

Fig. [3] shows the complex (and real) band structure 
in the [001] (fc z ) direction for two different values of feii: 

fen = and fc|| = (ko,0); the latter corresponds to one 
of the aforementioned "satellite" conductance hot spots. 
Arrows A and B indicate the pocket positions, and the ef- 
fective mass anisotropy is evident. The indicated pocket 
in the left panel (arrow A) has its long axis along the 
k z direction. The complex band (in orange, arrow A) 
departing from the edge of this pocket has a small curva- 
ture (as it is along the long axis of the pocket), leading 
quickly to high values of k for Ep < E c . On the other 
hand, the pocket indicated in the right panel (arrow B) 
has its long axis oriented in the k x direction. Here, the 
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FIG. 3: (color online) Real and complex band structure 
of Si along the k z direction at ku — (left) and fey = 
(0.85,0) (2ir/a) (right). Black lines stand for real bands as 
E — E(k z ); orange (or grey) lines stand for complex bands 
as E = E(k z + in). The effective mass anisotropy at the six 
pockets of the conduction band edge, k X]VtZ m ±0.85 (2n/a) 
evidently results in a quite different curvature of the complex 
bands departing from these (indicated by arrows). Therefore, 
when Ef is close to the conduction band edge, the tunnel- 
ing is dominated by the contributions at the four equivalent 
points fc|| w (±0.85,0) (2n/a), fey w (0, ±0.85) (2%/a). The 
symmetry of the bands is also indicated (Ai and As). (More 
complex bands, irrelevant for the discussion, have been omit- 
ted.) 



complex band along k z inherits a large curvature from 
the small axis of the pocket, and the value of k remains 
relatively small for Ep < E c . 

Consequently, when Ep is in the gap but close to the 
conduction band edge, the tunneling current through the 
satellite positions dominates, while the tunneling current 
through the Brillouin zone center is small. However, if 
Ep is lowered closer to the valence band edge, a different 
contribution from the Brillouin zone center becomes more 
important, marked by arrow C in Fig. [3] This comes 
about via the complex band of Ai symmetry arising from 
the valence band maximum. 

The picture becomes more complete if the full (001) 
surface Brillouin zone is scanned for the lowest branch 
of n(kf,Ep) at different values of Ep. Fig. [5] shows 
two such "complex Fermi surfaces": (A) one for Ep = 
E v ± 0.725 eV (close to E c ), and (B) one for E F = 
E v ± 0.325 eV (just a little lower than the middle of the 
gap). Brighter colored regions correspond to lower values 
of k. In the first case, as Ep is close to E c , the parts of 
the complex Fermi surface at the satellite hot spots show 
the lowest k. As we shall see in the next section, conduc- 
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FIG. 4: (color online) Schematic representation of the six 
equivalent Si conduction band "pockets". When these are 
projected on the two-dimensional (001) surface Brillouin zone, 
five hot spots for conductance in the [001] direction arise: one 
central (at feii = 0) and four equivalent satellites. 



tance calculations also show that these regions dominate 
the tunneling current for an analogous choice of Ep. In 
the second case shown in Fig. O the minimum of k is 
found at the Brillouin zone center. Then the tunneling 
contribution at fen = dominates if the Schottky barrier 
is thick. 

We close this section with the conclusion that the 
indirect gap of Si lends features to the complex band 
structure which are qualitatively different than direct- 
gap compounds, as GaAs, ZnSe, or MgO. Contrary to 
all these direct-gap materials, where the complex band 
of Ai symmetry at the Brillouin zone center gives the 
dominant contribution irrespective of the exact position 
of Ep, in Si the positioning of Ep can make a stark differ- 
ence. This effect can have consequences for all physical 
properties which depend on the complex band structure. 
Particularly in spintronics applications it can affect spin 
injection, tunneling magnetoresistance, but also ground- 
state properties such as inter layer exchange coupling. 24 



V. TOTAL AND SPIN-DEPENDENT 
CONDUCTANCE. CURRENT POLARIZATION 

We now proceed to the presentation and discussion of 
the ab initio results on the conductance and current po- 
larization. We start by commenting on the effect of the 
underestimation of the gap in density functional theory. 
The conclusions of Sec. lIVI are qualitatively still valid, but 
now, within the barrier, the possible choice of E c — Ep is 
more limited (otherwise Ep will enter the valence band). 
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FIG. 5: (color online) Complex band structure (decay param- 
eter ft(fc||; Ef)) in the band gap of Si in the [001] direction for 
E F = E v + 0.725 eV (A) and at E F = E v + 0.325 eV (B). In 
the former case, Ef is close to the conduction band edge, and 
the four equivalent contributions at fen ~ (±0.85,0) (2n/a) 

and fen ~ (0, ±0.85) (2n/a) dominate the tunneling; in the 
latter case, Ef is located a little lower than tne middle of 
the gap, and the decay parameter at fcy = is lowest. Only 
the lowest branch of K(k^) is shown. Brighter shaded regions 
correspond to lower k and more efficient tunneling. 



This results in an overestimation of the relative contribu- 
tion of the satellites to the current. As the calculations 
show, the satellites dominate the current even when Ep 
is at the mid-gap position; hence, calculated results at 
mid-gap are expected to be closer to a realistic situation 
of shallow tunneling close to E c . An overestimated prox- 
imity of Ep to the the valence band is necessary so for 
the Brillouin zone center contribution to prevail. More- 
over, the decay parameters and exponential falloff of the 
conductance with barrier thickness are underestimated. 



A. Total conductance 

Fig. [S] shows the calculated conductance per two- 
dimensional interface unit cell as a function of the barrier 
thickness. Results for two different positions of Ep are 
presented: mid-gap tunneling (Ep = E v + 190 meV) 
and shallow tunneling close to the valence band edge 
(Ep = E v ± 28 meV). The partial contributions of the 
region around the Brillouin zone center and the satellites 
are also shown. An exponential decay with thickness is 
characteristic of tunneling, and the slope of the curve in a 
logarithmic scale is proportional to the decay parameter. 
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FIG. 6: (color online) Conductance per two-dimensional unit 
cell as a function of the Schottky barrier thickness for (A) mid- 
gap tunneling (Etunn = 190 meV) and (B) shallow tunneling 
close to the valence band (Etunn — 28 meV). The contribu- 
tions of the region around the SBZ center (integrated around 
the f-point) and the satellites are shown separately. 



In the case of mid-gap tunneling (Fig. EJA_) it is clear 
that the satellite contribution dominates at all thick- 
nesses, being practically equal to the total conductance, 
while the central contribution to the total current is in- 
significant. While the slope of the satellite part of the 
conductance is constant, reflecting a single decay parame- 
ter K sat , the central contribution gradually changes slope; 
after a thickness of about 20 layers, only the Ai contri- 
bution remains around fc|| = 0, with a decay parameter 
Ki ~ k sa t- 

The behavior in the case of shallow tunneling, closer 
to the valence band edge (Fig. 03), is analogous. Here, 
however, the low positioning of Ep evidently results in 
Ki < K sat (see also Fig. Thus after a thickness of 
35 layers the central contribution prevails and dominates 
the tunneling current. 

In both cases we also observe that, for small thicknesses 
(5 layers), the satellite contribution is much stronger than 
the central one. Apparently there is a much stronger 
coupling between the Fe and Si states at the satellite 
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positions in the Brillouin zone. At this point we have no 
intuitive explanation for this effect. 



B. Spin dependent conductance and current 
polarization 

The spin polarization of the current is defined as 



P = 



h + h 



(i) 



where I-\ and Jj are the current of majority-spin and 
minority-spin carriers, respectively. The difference be- 
tween J| and arises mainly from the spin-dependent 
scattering at the Fc/Si interface (which gives rise to spin- 
dependent tunneling), due to the difference in coupling 
of Fe wavefunctions of different spin to Si wavefunctions 
at the interface. 

The calculated current polarization as a function of 
barrier thickness is shown in Fig. [7] Again, two cases 
are presented, corresponding to (A) mid-gap tunneling 
(Ep = E v + 190 meV within the barrier) and (B) shal- 
low tunneling close to the valence band edge (Ep = 
E v + 28 meV within the barrier), as was the case in 
Fig. [51 These were chosen as representative of differ- 
ent physical situations, where the dominant contribution 
to the current stems from different parts of the Brillouin 
zone. In both cases the band alignment in the injec- 
tion region (i.e., after the barrier) was chosen such that 
5 := Ep — E c = 28 meV. Apart from the total cur- 
rent polarization, two contributions of special interest 
are shown: one at exactly the f -point (P r ) and one in- 
tegrated around the r-point (P cn tr)- In terms of the k\\- 
resolved current, I-\(k\\) and Ii(kn), these contributions 
arc defined as 



Pr 



J T (fc|| =0)-/ i (fc|| =0) 
J T (fe|l =0) + / i (fc|| =0) 



(2) 



and 



Pcn.tr — 





i T (fcn) 




/ d2 h 







(3) 



where the latter integration takes place in the central 
part of the Brillouin zone where the current is non-zero 
(the central hot spot of Fig. [4j) . 

We first discuss Fig.[7]A Here the polarization is domi- 
nated by the satellite contributions for all thicknesses, as 
the current at the SBZ center is negligible (cf. Fig.[6K). 
The polarization is rather insensitive to the barrier thick- 
ness, being around P = 60%. In Fig.[7C, the fc||-resolved 
conductance is shown for both spin directions in the full 
SBZ for a barrier thickness of = 40 monolayers. 
Evidently the SBZ center has a negligible contribution, 
while the conductance at the satellite positions is higher 



for majority spin than for minority spin by approximately 
a factor of two. However, the contribution at exactly the 
r-point shows an interesting behavior, almost reaching 
the ideal P r = 100%; the integrated value around T, 
P = Pcntr, is lower. Before analyzing this we discuss 
Fig. [33. Here, for small barrier thicknesses, the current 
is dominated by the satellites, which yield a polarization 
of about 50%. But for larger thicknesses, only the central 
contribution to the current is of significance (cf. Fig. [6)3). 
Already at a thickness of 50 monolayers, P = P cn tr- And 
again we see that, exactly at the SBZ center, P r k, 100%. 

In order to elucidate the situation at the SBZ center we 
show in Fig. 03 the fc|| -resolved conductance focused in 
the region around the Y point, for the case of shallow tun- 
neling (Ep = E v + 28 meV in the barrier) and a barrier 
thickness of 60 monolayers. The majority-spin conduc- 
tance shows a single-peak structure with a maximum at 
the r point. On the other hand, the minority-spin con- 
ductance as a function of fc|| has a double-peak structure 
around k\\ = 0, with a pronounced dip at exactly fc|| = 
which is responsible for the value Pp « 100%. 

Obviously the r-point enjoys special properties. This 
"F-point rule" has been observed and explained in previ- 
ous work o 18 i 19 ' 25 i 26 on spin injection from Fe into direct 
gap semiconductors in the zincblcnde structure (GaAs, 
ZnSe, InAs). The reason is traced back to the symmetry 
of the wavefunctions at Ep, at fc|| = 0. The semiconduc- 
tor wavefunctions at fcj| = have Ai symmetry, which for 
Fe is present among the majority spin wavefunctions but 
absent among the minority spin wavefunctions at Ep, at 
least in the [001] direction. This symmetry mismatch of 
the Fe minority spin to the semiconductor wavefunctions 
results in almost total reflection, so that the current is 
almost 100% polarized. Departing from fc|| = 0, the Si 
bands acquire a mixed character, such that the Fe A5 
states (of d xz and d yz character), coupling to the Si A5 
complex band in the barrier, can also tunnel into the con- 
duction band after the barrier. Then the minority-spin 
transmission rises, as is shown in Fig.[7jDi21 Note that no 
such special point appears in the satellite hot spots al- 
though their centers lie on the high-symmetry directions 
r — M (along the cubic x and y axes). 

This behavior close to f is typical also for smaller 
thicknesses and for different tunneling energies, therefore 
the polarization P r reaches high values also in Fig. [7JA 
We infer that the integrated spin polarization can be in- 
creased if two requirements are fulfilled: (i) the central 
hot spot must be as small as possible; and (ii) the satel- 
lite contributions must be made negligible. Both can, in 
principle, be fulfilled, as we now discuss. 

Concerning requirement (i) , the radius k max of the hot 
spots depends on the injection energy S = Ep — E c in the 
injection region as S ~ But S is adjustable, e.g. by 

tuning the gate voltage of the doping concentration. In 
particular for Si, due to the high effective mass the value 
S = 28 meV used in the calculations hitherto is rather 
high and was chosen in order to reveal the structure 
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FIG. 7: (color online) (A and B) Current polarization as a function of the barrier thickness for an injection energy of 6 — 28 meV 
and for a tunneling energy Etunn at mid-gap (A) and at 28 meV (B). Shown are the total polarization (see Eq.[T]), polarization 
at F (Eq. [2j and the integrated polarization at the central hot spot (Eq. [3}- (C): Conductance as a function of k« in the SBZ 
for a barrier thickness of 40 layers in the case of Etunn at mid-gap. The SBZ center contribution is insignificant next to the 
satellites' contribution. (D): Similar as in (C), but for a barrier thickness of 60 layers in the case of Etunn = 28 meV. Here the 
0.1 x 0.1 central part of the SBZ is focused on, since the satellites are negligible. A double-peak structure, with a dip at T, is 
evident for minority-spin. 



of the fcii-resolved conductance, as already commented 
in Sec. [Til By choosing a smaller 5, the SBZ center is 
approached more and more, and the integrated polar- 
ization P cn tr rises. This is demonstrated in Fig. [8] At 
5 = 1.4 meV, P cn tr is already over 90%, while in the limit 



we obtain P a 



Pf w 100%. Interestingly this 



results in a "zero-gate anomaly" (if 5 is considered to be 
a gate voltage), demonstrated in the inset of Fig. [5] where 
Pcntr is shown as a function of J at a barrier thickness 
of 60 layers. Evidently the polarization drops abruptly 
with increasing S. Note that controlled injection at about 
S 1 meV requires low temperatures to avoid thermal 
broadening, since 1 meV corresponds to 11.6 K. In this 
respect, the T-point rule is easier satisfied in direct-gap 
semiconductors, where due to the low effective mass the 
fc|| -resolved current is confined to a tiny region around T 
also at higher 5. 

Although Fig. [5] shows results for shallow tunneling at 
Etunn — 28 meV, the behavior of P cn tr is entirely anal- 
ogous at higher E tunn , e.g., in the case of mid-gap tun- 
neling. Qualitatively, what changes is only the relative 
importance of the satellite contributions to the central 
contribution. 

This brings us to the discussion of requirement (ii). 
Is it possible to cut off the satellites? We already saw 
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FIG. 8: (color online) Polarization around the f point, Pcntr 
[see Eq. ©J, for small values of the injection energy <5. Inset: 
Zero-gate "anomaly" of the polarization for a 60-layer barrier. 



that this happens as the Fermi level approaches the va- 
lence band in the barrier region, if the barrier is thick 
enough. But there is also another possibility, namely by 
tctragonally straining Si. In tctragonally strained Si the 
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degeneracy of the six conduction pockets is lifted. In the 
case that c/a < 1, the two pockets along the c axis (i.e., z 
axis along the epitaxial direction in the presented geom- 
etry) are lowered in energy compared to the four pockets 
along the x and y axes^ Even a moderate strain of 0.5% 
(c/a = 0.995) was calculated to lift the degeneracy by ap- 
proximately 30 meV. As long as the injection energy S is 
kept under this limit, only the central conductance hot 
spot will be populated in the injection region, while the 
satellites will be cut off. 



VI. SUMMARY AND CONCLUDING 
REMARKS. OUTLOOK 

Electrical injection from Fe into Si(001) through a 
Schottky barrier has been shown to be theoretically pos- 
sible. A detailed discussion of the complex band struc- 
ture of Si in the [001] direction has revealed qualitative 
differences in the tunneling process compared to direct- 
gap semiconductors. As a result of the complex band 
structure, the current and polarization contributions at 
different conduction pockets have been shown to vary 
very strongly depending on the Schottky barrier thick- 
ness and the position of the Fermi level in the barrier. 
Depending on these parameters, the injection efficiency 
has been found to range between 50% and 100%. 

In the calculations a particular junction setup was as- 
sumed, including a number of approximations or idealiza- 
tions. One approximation lies in the shape of the Schot- 
tky barrier. It was assumed that the transition from the 
Schottky region to the injection region is abrupt. How- 
ever, but for the case of very precise interface engineering 
(e.g. with an appropriate doping profile), the transition 
to the injection region is more gradual. This would re- 
sult in position-dependent decay parameters, with the 
central part of the SBZ providing better tunneling close 
to the interface (where the middle of the gap should be 
at Ep) and the satellite positions being more efficient 
close to the injection region (where E c is lowered towards 
Ep). Thus, overall, either the SBZ center or the satellites 
would dominate the tunneling current, depending on the 
exact shape of the barrier. 



An idealization was that of an atomically flat Fc/Si 
interface, with the in-plane lattice structure unaltered. 
While possible in principle, in practice it can prove hard 
to achieve. If the two-dimensional periodicity is vio- 
lated at the interface, then the most severe consequence 
(as regards the results presented in the present work) 
is expected to be the absence of excellent spin filter- 
ing at the T-point. As was mentioned earlier, the ex- 
treme current polarization stems from the symmetry mis- 
match of the Si and minority-spin Fe wavefunctions at 
fcii = 0. Such symmetry arguments do not hold any 
more in the absence of perfect interface epitaxy. It 
has been shown^ e.g., in Fc/InAs(001) spin injection 
(where the same principle holds), that increasing inter- 
face disorder leads to a decrease of current polarization. 
However, a few perfectly epitaxial Fe layers should be 
enough for a symmetry-induced polarization. This has 
been found in an analogous case of symmetry-induced 
polarization in Fe/MgO/Fe(001) tunnel junctions.— Fur- 
thermore, it should be noted that, since an MgO barrier 
is known to be selective of the Ai states at fc|| = 0, it 
can also be used to increase the efficiency in spin injec- 
tion experiments. Efficient spin injection has been found 
in FcCo/MgO/GaAs and FcPt/MgO/GaAs junctions^ 
while work in this direction has been reported also for 
FeCoB/MgO/Si junctions.- 1 ^. 

A reduction of efficiency can also be caused by the for- 
mation of iron silicide at the interface, which can some- 
times lead to non-collinear magnetic ordering. This can 
be avoided by inserting a nonmagnetic metal between Fe 
and Si (as was done, e.g., in Refs. [9IIT0I): calculations on 
such junctions will be the object of future work. 
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